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Distortions of the \/3 x \/3 Sn/Ge(lll) and Sn/Si(lll) surfaces are shown to reflect a disproportion- 
ation of an integer pseudocharge, Q, related to the surface band occupancy. A novel understanding 
of the (3 X 3)-lU ("1 up, 2 down") and 2U ("2 up, 1 down") distortions of Sn/Ge(lll) is obtained 
by a theoretical study of the phase diagram under strain. Positive strain keeps the unstrained value 
(5 = 3 but removes distorsions. Negative strain attracts pseudocharge from the valence band causing 
first a (3 X 3)-2U distortion (Q = 4) on both Sn/Ge and Sn/Si, and eventually a (V3 x ^/3)-3U ("aU 
up") state with Q = &. The possibility of a fluctuating phase in unstrained Sn/Si(lll) is discussed. 



Disproportionation, namely the spontaneous differen- 
tiation of initially equivalent sites into states of unlike va- 
lence or charge, is well known in solid state chemistry, for 
example in organic charge-transfer saltall and also in tran.- 
sition metal and mixed valence rare earth compoundsa. 
Common features are an enlarged unit cell, due to or- 
dering of the disproportionated sites, and the presence of 
a well defined quantum number (e.g., the valence) that 
is integer, even when the actual charge inbalance is gen- 
erally not integer. Here we propose that a dispropor- 
tionation, driven in this case by a different physics, but 
equally characterized by integer valencies, is at the ori- 
gin of the well known but only ip.completely understood 
two-dimensional phase transitiona observed in 1/3 mono- 
layer Sn-covered Ge(lll) and Si(lll) surfaces. Moreover 
we present theoretical suggestions that the disproportion- 
ated state could be changed, or even removed altogether, 
by external strain. 

The 1/3 monolayer Sn-covered Ge(lll) and Si(lll) 
surfaces have received considerable attention since the 
discovery of a transition from a room-temperature phase 
with \/3 X V3^30° periodicity to a low-temperature 3; 
phase, first seen in Pb/Ge(lll), then in Sn/Ge(lll)E , 
and unsuccessfully searched for in Sn/Si(lll)B. The main 
reason of interest lies in the electronic driving force be- 
hind the transition, naniely the presence of a half filled 
surface state, suggestive^ of a two-dimensional surface 
charge-density wave (SCDW) mechanismEl. Recent X-ray 
data show that the transition is accompanied by static 
freezing and long-range ordering of a 3 x 3 periodic struc- 
tural distortion, consisting mainly of a "bucklingjl^ x)f Sn 
adatoms by ~ 0.35 A along the surface normalD'cl. In 
the distortion, adatoms are generally believed to follow a 
"1 up, 2 down" (lU) pattern, but there is evidence that 
under some as yet unclear circumstances a non-negligible 
portion of the surface distorts with th&,p»posite phase, 
namely a "2 up, 1 down" (2U) patternBllj. The precise 



mechanism underlying the phase transition has led to de- 
bate. While the driving role of the partly filled surface 
states is unquestioned, the 2D Fermi surface area and 
nesting properties of Sn/Ge(lll)j_d£L.not fit a 3 x 3 pe- 
riodicity as expected for SCDWsEJ'Ej. Due to the large 
electron-electron repulsion, a Mott insulator or a spin- 
density wave stateEj with negligible distortion could pre- 
vail, as in the isoelectronic case of SiC(0001)lia. The 
large observed distortion and the order-disorder nature 
of the phase transitiorO however rule out both the weak- 
coupling SCDW and the Mott insulator. First-princ 
electronic structure calculations, initvaJUi, negativel 
later correctly favored the lU distortiorJla^tj, with siirfiice 
electronic bands in agreement with experimentu'O LI . 
The calculated energy gain was about 10 meV/adatomlla, 
mostly of Jahn- Teller nature. Multi-band and bond- 
density wave aspects were also pointed out as crucial fac- 
tors that|-sjtabilize the lU distortion relative to the Mott 
insulatoJla. Surprisingly, these-|Calculations showed no 
energy gain for a 2U distortionlla, contrary to the kna|iai 
readiness of 2U to show up experimentally in Sn/GeoE2l. 
The order-disorder character of the phase transition, well 
confirmed by simulationsEil, is in agreement with the per- 
sistence of a lU-like Sn Ad core level splittingj-Ln nomi- 
nally undistorted Sn/Ge at high temperaturcEj. How- 
ever a puzzling 2U-like uersistent core splitting in appar- 
ently undistorted Sn/SitZfE3 remains unexplained. Fi- 
nally, substitutional defects play an important role. A 
Ge (Si) adatom replacing Sn acts as an electron donor 
in the Sn/Ge (Sn/Si) susface, triggering a local distor- 
tion of roughly 2U typecl. However the defect is also 
likely to exert a negative local strain, effectively squeez- 
ing the neighboring Sn adatoms outwards relative to the 
substrate, so that strain and charge donation apparently 
conspire. One is eventually led to wonder about the gen- 
eral role of strain, and about whether some of these open 
questions could be addressed by studying primarily the 
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effect of strain on these surfaces. We report here the sur- 
prising results of a theoretical study of the Sn/Ge(lll) 
and Sn/Si(lll) systems under uniform strain. Besides 
predicting new strain-driven phase transitions, we find 
that these surface distortions are in the nature of dispro- 
portionation phenomena. 

We carried out electracLic structure calculations within 
the gradient correctedEEl local density approximation 
(GC-LDA) where, in the absence of strain, a stable 
(3 X 3)|iLUiistortion was established for the Sn/Ge(lll) 
surfacdi3'E3j— I whereas Sn/Si(lll) is still marginally 
undistortcdEil. We used slab geometries consisting of 
three Ge (Si) bilayers, a planar 3x3 cell with 18 
atoms/bilayer, 3 Sn adatoms per cell in the T4 position, 
and 9 hydrogen atoms saturating the bottom surfae© 
of the slab. We employed the PWSCF parallel codeE2l 
with a well checked plane-wave cutoff of 12 RyEl, and It- 
summation schemes based on 6 Monkhorst-Pack k-points 
in the irreducible zone, extended when necessary to 16 
and 36 k-points. We strained the slab cell isotropically 
in the surface plane, by increasing or decreasing its (x, y) 
sides by up to ±5%, full relaxation allowed along z, the 
surface normal, so that inside the slab the three dimen- 
sional strain tensor was effectively uniaxial. For each 
strain value, the stable surface geometry was sought for 
through a series of relaxation runs performed by con- 
straining the value of the adatom buckling Az, (the z- 
difference between "up" and "down" adatoms) and relax- 
ing all other atomic coordinates, except for the lowest 
Ge (Si) monolayer, held fixed and bulk-like (or strained), 
along with the saturating H layer. 

Figure |l] summarizes the zero-temperature phase dia- 
gram versus strain, the order parameter and energy gain 
being given for the 3x3 distorted state relative to the cor- 
responding fully relaxed \/3 x -v/S state at the same strain, 
for both Sn/Ge and Sn/Si. At zero strain, in agceemettt 
with experiment and with previous calculational3'Ej'ta, 
Sn/Ge is lU distorted, with Az ^ 0.34 A, whereas Sn/Si 
is OU (undistorted) . Positive strain (in-plane expansion) 
causes Sn/Ge to lose its distortion, through a lU-OU 
transition. By mechanically attracting the Sn adatoms 
closer to the Ge substrate, strain pushes the half filled 
dangling bond band upwards in energy towards the mid- 
dle of the gap. In line with our earlien finding that the 
3x3 distortion is a multi-band effectll3, in this restored 
single-band situation the surface reverts to the undis- 
torted V3x V3 state. The resulting narrow half-filled sur- 
face band metallic state should actually be prone to other 
instabilities. We did check with spin-polarized calcula- 
tions for a possible antiferromagnetic (Mott-Hubbard) 
state similar to that of isoelectronic SiC(OOOl); but up 
to 5% strain we did not find one, at least in Sn/Ge. 
Negative strain (in-plane compression) on the contrary 
pushes the adatoms upwards and that, beyond a criti- 
cal strain of ~ —3%, stabilizes a 2U distortion, with a 
1U-2U transition in Sn/Ge, and a 0U-2U transition in 
Sn/Si. Within our approximation, the 2U order param- 
eter Az at e = —4% is 0.49 A with an energy gain with 



respect to the OU configuration of 39.5 meV/adatom in 
Sn/Ge, and 0.4 A with an energy gain of 30 meV/adatom 
in Sn/Si. In both cases the energy gain is suggestive 
of a 2U-disordered transition at relatively large critical 
temperature. At still larger negative strains ~ —5%, all 
adatoms are finally forced up, with a 2U-3U transition, 
and an optimal 3U state for both Sn/Ge, and Sn/Si. This 
state should be seen as uniformly distorted - same sym- 
metry as the "OU" undistorted ^/S x Vs state, but with 
quite different physics as we shall see. 

Closer examination of the distorted geometries reveals 
a striking feature: the Sn adatoms' height above the sub- 
strate is bimodal, approximately quantized for all strains 
and states. A given Sn adatom can be either U or D, 
but not in between. This is reminiscent of static mixed 
valence, or disproportionated states, where an individual 
ion site is endowed with a choice between two integer 
valence states, the overall proportion between them con- 
trolled by a chemical potential. The valence state can 
either ffuctuate dynamically with time, or be frozen in a 
static regular or irregular spatial pattern. In that case, all 
local properties, including the surrounding ion geometry, 
display a bimodal distribution. We found, by analysing 
the electronic structure, that in our case the role of ef- 
fective valence is played by a surface state occupancy q, 
which is 2 in all U adatoms, 1 in all adatoms of the undis- 
torted ("OU") phase, 0.5 in D adatoms of the lU phase 
of Sn/Gc, and zero in the D adatoms of 2U phases. This 
conclusion was established with the help of a separate 
Sn/Ge calculation, where a single Sn adatom was con- 
fined in the 3x3 slab cell. The six remaining first-layer 
Ge atoms were saturated by six hydrogens, effectively de- 
coupling the Sn adatom from its replicas. When a vari- 
able electron number was introduced with a net electron 
excess between —1 and -1-1 (charge neutrality obtained 
through a matching uniform positive background), it was 
found that the population change was restricted to the 
Sn dangling bond surface state, its occupancy q corre- 
spondingly varying from zero to 2. For each q value, 
the total energy was minimized as a function of adatom 
height z over the substrate, and the optimal 20(9) was 
determined. The calculations were repeated for various 
strains and the results reported in Fig. || (dashed fines). 
The comparison of the calculated zo((7)'s with the indi- 
vidual adatom heigths previuosly calculated in the OU, 
lU, 2U and 3U stable phases of Sn/Ge is striking and 
provides a strong evidence of quantization of the individ- 
ual (7i's at four and only four heigth values, corresponding 
within a few percent error to the occupancies 2, 1, 1/2, 
and indicated above. These quantized charges qi are ac- 
tually pseudocharges, not to be confused with real, phys- 
ical electron charges on the adatoms. An approximate 
MuUiken analysis of valence chargeaHa showed in fact that 
the actual charge redistribution among adatoms is always 
less than 0.1 e. A negligible redistribution of true, total 
charge is also in good accord with the observed Ad Sn core 
level splitting, consisting of one deep and one shallow lev- 
els in a 1:2 intensity proportion, attributable to U and 
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D adatoms respectively. Had the true charge inbalance 
even only resembled that of the pseudocharges, 2 and 0.5 
for U and D adatoms, respectively, the sign of the split- 
ting should have been the ppposite, and its magnitude 
much larger. We calculatcdnJ the magnitude of the core 
level splitting between U and D adatoms in Sn/Gc-lU in 
the final-statec^ approximation, and obtained 0.20 eV, 
in reasonable agreement with 0.32 eV, observed both in 
(3 X 3)-lU and-fotlbe high temperature -s/S x V3 phase 
of Sn/Ge(lll)lljll3'Ea. The split core level intensity ratio 
should be instead 2:1 for a 2U distorted surface, the U 
and D adatoms being now in that proportion, opposite 
to the 1:2 of the lU case. Repeating our core level calcu- 
lation for the -3% strained Sn/Si(lll) (3 x 3)-2U surface 
we obtained a Sn Ad splitting of —0.20 eV, in qualitative 
agreement with ~ —0.4 eV observed on the nominally 
undistorted Sn/Si(lll) V3 x V3 surface. We shall return 
to this point below. 

We are now in a position to address directly the physics 
behind the pseudocharge and its disproportionation, with 
and without strain, in connection with the electronic 
structure. In the OU -x/S x phases, the surface state 
band is isolated inside the gap, each adatom contribut- 
ing one electron to it, whence qi = \ and Q = Y^qi = 3. 
This is an undistorted metal state, at least within GC- 
LDA. In the (3 x 3)-lU phase, the surface band hybridizes 
much more with the valence bands, especially with the 
Ge-Ge bonding state directly underneath the adatomtj. 
There is still on average one electron/adatom; however 
the total pseudocharge (5 = 3 per cell disproportionates, 
q — 2 electrons now filling the U adatom's dangling bond 
(Figure ^) , the remaining one shared between the other 
two D adatoms, each with q — 0.5. The onset of the 
(3 X 3)-lU ground state in Sn/Ge is thus favored by the 
small Ge bandgap, and by the consequent energy prox- 
imity of the surface band to bulk valence bands. That is 
consistent with the absence of distortion in Sn/Si(lll), 
and in SiG(OOOl), where the surface band is more iso- 
lated in the wider bandgap. In the strain-induced 2U 
phase, both U adatoms have q = 2, while D has q — 0. 
The total pseudocharge Q is now 4: one extra electron 
has been drawn from bulk valence bands into the sur- 
face band. For strain values in the 2U range, the surface 
bands (Figure ^) drop in fact enough in energy to over- 
lap the bulk valence bands, now slightly above Ep- Fur- 
thermore, because the difference qu—qo is greater than in 
the lU phase (2 compared to 1.5), the optimal distortion 
is larger (0.49 A compared to 0.34 A in Sn/Ge). The 2U 
distortion amplitudes and energy gains are correspond- 
ingly smaller in the stiffcr, wide-gapped Sn/Si. Finally, 
in the (\/3 x -\/3)-3U phase, all adatoms are equivalent 
with q — 2. The surface band lies entirely below Ep, 
and 3 electrons are transferred from bulk to surface, so 
that (5 = 6. In this sense, the 3U state is fundamen- 
tally different from the OU state, despite its sharing the 
same symmetry. In all these distorted states, the sur- 
face state band breaks up into narrow adatom-localized 



levels. Their filling and emptying leads to a Jahn- Teller- 
like energy gain, and that explains the quantization of 
the pseudocharge q. Negative strain acts as a chemical 
potential, and disproportionation arises out of the neces- 
sity to distribute an increasing electron number through 
integer surface state occupancies. In chemical terms, U 
adatoms could perhapsJpe compared with "tetravalent" 
Sn, D to "divalent" SnEZl. 

These results and predictions should be amenable to 
experimental test, particularly if the strains required to 
provoke the transitions could realistically be introduced, 
e.g., by heteroepitaxy. Even if our calculated critical 
strain values are probably affected by systematic errors of 
order 1-2%, connected with lattice parameter errors and 
with imperfect handling of incipient strong correlations 
by the density functional calculation, the strain-induced 
sequence of transitions seems inescapable. Among the 
phenomena that ape already observed, the 2U patches 
reported on Sn/GeO seem now very likely to be induced 
by some local strain. The puzzling 2U core level signa- 
ture of nominally unstrained Sn/Si is likely to be due to 
a dynamical, fluctuating 2U state. Assuming in fact a 
0U-2U zero-temperature, classical transition to occur in 
reality at some small positive strain (instead of a small 
negative strain as in the GC-LDA calculation), the long 
range order could simply be destroyed by thermal fltictu- 
ations, still leaving a local 2U distorted structurecS. Of 
additional interest is the possibility of quantum mechan- 
ical fluctuating states, analogous to dynamical mixed 
valencea. Glose to the critical strain values separating 
two classical distortion states, such as OU-IU, 0U-2U, 
1U-2U, 2U-3U, tunneling across low barriers could lead to 
quantum melting of the classical long-range order. This 
kind of polaronic state could entail the exotic but real 
possibility of surface state superconductivity when, as the 
case for 0U-2U, the pseudocharge fluctuation takes place 
in units of two, thus realizing a bipolaron state. The 
large distortions involved suggest however that this novel 
strong-coupling surface polaronic metals and bipolaronic 
superconductors might be attainable only at very low 
temperatures. 
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In this kind of dynamical state, as in mixed valence, Sn 
adatom properties in Sn/Si should appear bimodal for fast 
probes, and unimodal for slow probes. Preliminary data 
indicate a single Sn site in photoelectron diffraction at 100 
K (S. Modesti, private communication). 
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FIG. 1. Energy gained by the (3x 3)-lU and 2U distortions 

against the relaxed undistortod a/S x y/3 phase, as a function 

of strain, (a) Sn/Ge(lll); (b) Sn/Si(lll). 
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FIG. 3. Electron bands of Sn/Ge(lll) 3x3. (a) lU, zero 
strain; (b) 2U, strain —4%. 
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FIG. 2. Adatom height z as a function of strain in 

Sn/Ge(lll). Open squares: "down" adatoms; filled squares: 

"up" adatoms; dashed lines: optimal height of an isolated 

charged adatom (see text). 
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